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Abstract 

We study the asymptotic behavior and the asymptotic stabihty of the two-dimensional Euler equations 
and of the two-dimensional linearized Euler equations close to parallel flows. We focus on flows with spectrally 
stable profiles U (y) and with stationary streamlines y = yo (such that U'{yo) = 0), a case that has not 
been studied previously. We describe a new dynamical phenomenon: the depletion of the vorticity at 
the stationary streamlines. An unexpected consequence, is that the velocity decays for large times with 
power laws, similarly to what happens in the case of the Orr mechanism for base flows without stationary 
streamlines. The asymptotic behaviors of velocity and the asymptotic profiles of vorticity are theoretically 
predicted and compared with direct numerical simulations. We argue on the asymptotic stability of this 
ensemble of flow profiles even in the absence of any dissipative mechanisms. 



o _ 

1 Introduction 



CsJ ■ The flow of a perfect fluid is described by the Euler equations, one of the oldest equations of mathematical 
^ I physics [28J. Four century after their discovery by Euler, these equations still propose big challenges both 
to mathematician and physicists |28| . Two-dimensional flows and the two-dimensional Euler equations are 
mathematically much simpler than their three dimensional counterparts, but still present some very interest- 
ing unsolved problems. One of the main phenomena for two-dimensional flows is the self-organization into 
coherent structures dH [Ml 1121 ISH 1111 112] : monopoles, dipoles, and parallel flows. Such large-scale structures 
are analogous to geophysical cyclones, anticyclones, and jets in the ocean and atmospheres. This analogy, 
understood thanks to the theoretical strong similarities between the 2D Euler equations on one hand and the 
Quasi-Geostrophic or the Shallow Water models on the other hand, is one of the main motivations for the study 
of the 2D Euler equations. The 2D Euler equations also describe experimental flows: the transverse dynamics 
of electron plasma columns [57j , the dynamics of fluids when three-dimensional motion is prevented by a strong 
. ^ transverse field (rotation, a transverse magnetic field in a liquid metal, etc. |64| ) or the dynamics of fluids in 
■ very thin geometries [51]. 

Because large-scale coherent flows appear spontaneously in two dimensional turbulence, their stability is a 
crucial problem. Moreover, the study of the dynamical mechanism that describes the relaxation towards these 
stable flows is essential. In this paper, we consider stable parallel base flows vo(x,?/) = U{y)ex, which are 
dynamical equilibria of the 2D-Euler equations. We prove that the velocity of these flows is asymptotically 
stable^ that is that all the solutions to the nonlinear Euler equations that start near Vq converge to some other 
parallel flows vq -1- 5U {y) e^^, close to vc0- Our analysis mostly relies on the linearization of the Euler equations 
close to the base flows, which we prove to actually describe also the nonlinear relaxation at leading order. More 



^This notion of convergence makes the notion of asymptotic stability stronger than the alternative Lyapounov stability, that 
only states that all solutions of the nonlinear equations that start near a steady point vo stay near vo forever. The asymptotic 
stability of the velocity refers to asymptotic stability in the kinetic energy norm. 
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precisely, we prove that the perturbation velocity decays algebraically for large times. As far as the linearized 
dynamics is concerned, an important improvement over the previous works is the understanding of the case when 
the flow has some stationary streamlines y — yo (or equivalently the velocity profile U {y) has some stationary 
points yo, U' (yo) — ojl, which has not been elucidated even qualitatively previously. 

Besides the stability and the asymptotic stability problem itself, the evolution operator for the linearized Eu- 
ler equations plays a very important role in different statistical approaches to turbulent flows [24l [211 [31] . Indeed, 
in a turbulent context, it is likely that the qualitative or quantitative properties of the fluctuations around such 
stable structures are related to the linearized dynamics. For instance, in quasi-linear approaches or second order 
closure of the evolution of the Euler equations, point vortex model [24J, or of the Navier-Stokes equation (either 
forced or unforced, and either deterministic or stochastic,) the linear operator appears naturally as an essential 
theoretical tool. Similarly, in the forced problems, in the linear regime, the response can be easily expressed in 
terms of the evolution operator for the linearized dynamics or of the resolvent operator for the dynamics. The 
behavior of such operators for large times is thus a very important issue that has many theoretical and dynamical 
consequences. In this work, we quantify very precisely the large time asymptotic behavior of the evolution op- 
erator for the linearized 2D-Euler equations, and discuss briefly some of the implications for the above problems. 

The stability of the large-scale coherent structures of two-dimensional flows is a very old and classical fleld 
of fluid mechanics. For instance Rayleigh [521, Kelvin [55], Orr [SO], Sommerfeld [HJ and many other famous 
scientists from the nineteenth and the beginning of the twentieth century have participated to the understanding 
of the linear theory for the 2D Euler equations close to parallel flows. Mathematicians gave also important 
contributions: Arnold's theorems [5] and some modern generalizations [551 ICTl [251 [T51 ITU] prove the Lyapounov 
stability of some of these flows. Even if this work deals only with the behavior of slightly perturbed stable flows, 
other equally interesting and important problems arise in the study of unstable or oscillatory flows. Recently 
many works have been devoted to the proof of the instability of some classes of flows, the characterization of 
the spectrum of the linearized equations, and some estimates on the stability and the instability of such flows ; 
see for instance 1501 1571 [5] and references therein. 

Our work is based on the linearized 2D Euler equations. The Rayleigh equation [53], which describes modes 
for the linear dynamics, has been a subject of mathematical and theoretical research since the beginning of 
the twentieth century |23| . and is still currently active. The main interest lies in the dynamical phenomena 
associated with the singularities at the critical layers (the singularities appearing when the frequency of the 
perturbation is equal to that of a closed streamline of the base flow). However, the modes of the Rayleigh 
equation do not describe the fully linearized dynamics, because the linear operator is non- normal [30J. Among 
the peculiarities of the linearized Euler equations, we stress the Orr mechanism |50J: the base flow shears 
the perturbation producing thinner and thinner fllaments ; then when the velocity or the streamfunction is 
computed, the effect of such fllaments being smoothed out, the perturbation velocity decays for large times. 
This mechanism is easily quantifled when the shear is linear, for the Euler equations [S^ or using for instance 
Kelvin waves [55] for viscous flows. Case [H] and Dikii [22] were the flrst to stress that in general, for inviscid 
flows, such a phenomenon is outside the scope of a modal description using the Rayleigh equation. When the 
shear is linear, using a Fourier-Laplace transform, the dynamics of the perturbation is properly described in the 
framework of an initial value problems. They concluded that the perturbation velocity decreases asymptotically 
with an algebraic law for large times. Other phenomena associated to the non-normality of the linear operator 
include possible transient growth [291 1551 [T]. inviscid damping (the counterpart of Landau damping in plasmas), 
axisymmetrization [56' '40] , and algebraic instabilities [49] . From a mathematical point of view, the singularities 
at the critical layers lead to the existence of a continuous spectrum for the linearized Euler equation. The 
analysis of the properties of this continuous spectrum explains most of these transient growth, inviscid damping, 
algebraic instabilities, and so on. 

From a theoretical point of view, one class of works used the Laplace transform tools [5] HI [57] , following 
the initial works of Case [19], Dikii [22] and the generalization to non- uniform shear by Briggs-Daugherty and 
Levy [TS]. Another class of studies, less general but very enlightening, used simple or particular base flows or 
special conditions for which explicit computations are possible |65| [501 HOI [IE [051 [HI [HI [IHl [IS] ■ 

In this paper, we are especially interested in the precise description of the large time asymptotic behavior 
of the 2D Euler and 2D linearized Euler equation close to parallel base flows. For the linearized dynamics of 
stable base flows, once the contribution of possible neutral modes has been subtracted, the asymptotic behavior 
is related to the continuous spectrum of the linearized operator. In the case of the base flow with a linear shear, 

^Please note that the profile U{y) has stationary points, but the two-dimensional base flow has no stationary points 
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U {y) = ay, the explicit computations by Case [19] showed that, due to the Orr mechanism, for large times the 
velocity perturbation decays algebraically : 
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with exponents a = 1 for the longitudinal component Vx of the velocity perturbation and /3 ~ 2 for the transverse 
one Vy. 

For more general base flows with strictly monotonic profiles U (y) (without stationary streamline), it is a 
common belief that the exponents a = I and (3 = 2 remain valid. This belief is based on the results of an ansatz 
for large time asymptotics [T3] (see also [H], appendix A). Some interesting comments about the temporal 
behavior of the streamfunction and velocity, in the case of localized initial perturbations (vorticity defects), can 
also be found in section 7. Even if we have not found any complete rigorous proofs, very precise classical 
arguments using the Laplace transform [SSI [H] conclude that the contribution of the continuous spectrum to 
the stream function perturbation -0 decays with ij; = O (j), in agreement with equation ([T]). However, these 

arguments do not generalize where the profile U{y) is not monotonic (flows with stationary streamlines). 

From Lundgren work ([H], appendix A), we see that the preceding algebraic decay for the velocity or the 
streamfunction may be related to the following asymptotic behavior for the perturbation vorticity: 



where k is the initial perturbation wave number. Indeed, computing the velocity from Lundgren's ansatz 
(equation ^) and assuming uniformity in the asymptotic expansion, we obtain oscillating integrals leading to 
algebraic decay for large times. The values for the exponents a = 1 and (3 = 2 are then related to the singular- 
ities of the Green function used in order to compute the velocity perturbation from the vorticity perturbation. 
This argument, assuming Lundgren's ansatz, suggests that the asymptotic behavior for the velocity should be 
different for velocity profiles U{y) with stationary points j/o {U' (j/o) = 0, base flow with stationary streamline). 
Actually, in such a case, the stationary phase asymptotics for oscillating integrals would generically give l/\/t 
contributions. It has then been noticed by several authors, that with such a law, Lundgren's ansatz would 
be self-consistent. [131 dHl [ID- Similar problems have also been noticed by Brown and Stewartson [T3], as 
their own asymptotic expansion clearly breaks down where U'{y) = (base flows with stationary streamlines). 
Besides Lundgren and Stewartson, many authors have insisted on the specificity of base flows with stationary 
streamlines (see for instance [TBI UM [55] ). 

In the past, there have been only a few studies considering base flows with stationary streamlines. In the 
case of the equations for 2D barotropic flows on a /3 plane (a direct generalization of the 2D Euler equations), 
Brunet and coauthors |16[ [H] have studied the dynamics close to a parabolic jet when the potential vorticity 
gradient exactly cancels the (3 effect. This case is similar to the linear shear case in the Euler equations studied 
by Kelvin, Orr. Case and others [SS] [SHI dHl E3 122 j in that the vorticity gradient exactly cancels out, which 
makes the linearized equation much simpler and amenable to a very interesting explicit analytic treatment. In 
the following, we will argue that, because of the cancelation of the vorticity (or potential vorticity) gradient, 
the dynamics of these cases is actually non generic, and that flows where the vorticity (or potential vorticity) 
gradient does not vanish behave differently. 

In the general case, the asymptotic behavior of the vorticity and velocity perturbations of flows with sta- 
tionary streamlines thus remains unstudied. In natural flows, however, jet velocity proflles are most of time 
not monotonic but have some extrema, i.e. the flow has stationary streamlines ; see for instance Jupiter, atmo- 
spheric, and ocean jets. Why have cases with stationary streamlines not been studied previously ? This may be 
partially on account of the wrong belief that base velocity profiles with stationary points should be unstable. It 
is true that many of the flows with extrema in their velocity profiles do not fulfill the classical Rayleigh-Fj0rtoft 
criteria [23J. However, these criteria are only sufficient conditions of stability. Moreover, as seen in natural flows 
and as shown bellow with several examples, many parallel flows with stationary streamlines and not fulfilling 
the classical Rayleigh-Fj0rtoft criteria are actually stable. Another reason for a lack of studies may also be the 
theoretical difficulty with Laplace tools in this case, related to the presence of stationary streamlines (merging 
of critical layers). Indeed, an essential tool for the Laplace transform is the analytic continuation of dynamical 
quantities, performed by avoiding the singularities associated to the critical layers, with the use of integration 
in the complex plane [23J. As will be discussed bellow, in the case of flows with stationary streamlines, in order 




(2) 
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to perform the analytic continuation, one would need to find a path in the complex plane passing at the same 
time above and bellow the singularity, which is clearly impossible. For this reason, it has often been stated that 
Laplace tools cannot be used when stationary streamlines are present. By contrast, we illustrate in this work 
that even if analytic continuation can not been performed in this way, Laplace tools are still very useful and 
lead to very interesting results. 

Li the following we consider the generic case of a parallel flow with any profile, either with or without 
stationary streamlines, improving by far the class of previously studied flows and overcoming the previously 
discussed difficulties. We also discuss possible generalizations to monopole vortices. We show how the Laplace 
transform is generalized to the case of base fiows with stationary streamlines. For instance, we show how 
the classical determination of the number of unstable modes, by using Nyquist's plots, remains valid in this 
case. From this general theoretical approach, we prove that the asymptotic vorticity field actually follows the 
Lundgren's ansatz ([5]), even in the case of a base fiow with stationary streamlines. Similarly the velocity field 
decreases also algebraically with the power laws ([T]), with a = 1 and /3 = This may seem paradoxical, 
after the discussion of the preceding paragraphs. Actually, the naturally expected contributions from the 

stationary phase asymptotics do not exist, unexpectedly. One reason is the non- uniformity of the asymptotic 
expansion in Lundgren's ansatz. Another more important reason is related to a new dynamical phenomena 
leading to the rapid decrease and cancelation of the vorticity perturbation exactly at the stationary streamline 
(see Fig. [T]), which partially erase the effect of the stationary phase. We call this phenomena vorticity depletion 
at the stationary streamlines. This is a non-local collective phenomena, due to the effect of the perturbation 
velocity on the background vorticity gradient. For this reason, this phenomena has not been observed in the 
previous studies involving stationary streamlines [161 [T¥l because these cases have an exactly zero vorticity 
(or potential vorticity) gradient. These last cases are thus non generic. 

In the following, we predict the vorticity depletion at the stationary streamlines using Laplace tools. It is 
thus a generic effect in any type of parallel flow with non-monotonic stable velocity profile. We also illustrate 
the results by direct numerical simulations in the case of Kolmogorov base flow U(jj) = cosy, for the 2D 
Euler equations with periodic boundary conditions. This vorticity depletion mechanism also impacts non linear 
turbulent flows when the perturbations are small enough to be governed by the linearized equations, as discussed 
in the conclusion. 

We establish the large time asymptotic contributions to the vorticity and to the velocity fields, by the con- 
tinuous spectrum, in the case of the linearized dynamics. We are then able to discuss the asymptotic stability 
of the velocity of parallel base flows, for the non-linear dynamics. The result is that parallel base flows which 
have no mode (neither stable nor unstable) are asymptotically stable: any small perturbation leads to a small 
deformation of the base flow, the perturbation velocity to this new base flow decaying algebraically. Stable 
Kolmogorov flows are examples of base flows without modes, illustrating the importance of this class of flows. 
Based on these results we also conclude that a quasi-linear approach predicts the asymptotic velocity profile. 
We note that all this is true only thanks to a non-trivial cancelation of leading order terms, already noted by 
[41| . The current work put the validity of a quasi-linear approach on a more rigorous ground, and proves that 
this it is also valid also for fiows with stationary streamlines. It also gives an efficient theory and numerical tool 
to predict the asymptotic flow. 

For the case of a circular vortex base flows, the work [7] shows that the far field velocity decays with ex- 
ponents different from the case of parallel flows. For circular vortex base flows, an example of non-monotonic 
angular velocity profile has also been studied [49] , based on a special explicit solution [62] . This example shows 
a very interesting algebraic instability with t^/^ growths. Even if we do not explicitly treat the case of stable 
circular vortices v (r) = [/ (r) eg in the present study, the generalization to this case of the present study could 
be taken following similar theoretical arguments. The mechanism of vorticity depletion at stationary streamlines 
or at the core of stable probably also exist for stable vortices, as discussed in the conclusion. 

In section[5]we introduce the 2D Euler equations and the linearized Euler equations. Section [2] describes the 
theory related to the linearized Euler equation. We discuss the main results related to the asymptotic behavior 
of the vorticity and velocity fields in section 12.21 The core of the proof relies on the results of the analysis on 
the limit of small e for the resolvent operator in section 12.31 

Section [3] discusses the asymptotic stability of parallel flows for the Euler equation. 

•^An exception is the velocity field close to the stationary streamline, where we have no theoretical prediction, but where we 
observe numerically that either a = 1 and /3 = 2ora = /3 = 3/2 depending on the symmetry of the perturbation. 
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Figure 1: Evolution of the vorticity perturbation lo{x, y,t) — uj {y, t) exp (ikx), close to a parallel flow vo(x, y) = 
U{y)ex with U (y) — cos (y). The figure shows the modulus of the perturbation \uj [y, t)\ as a function of time 
and y. One clearly sees that the vorticity perturbation rapidly converges to zero close to the points where the 
velocity profile U (y) has extrema (y = and tt) . This depletion of the perturbation vorticity at the stationary 
streamlines is a new generic self-consistent mechanism, understood mathematically as the regularization of the 
critical layer singularities at the edge of the continuous spectrum. 
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These results are illustrated in the case of doubly periodic boundary conditions, with the Kolmogorov flows 
as a base flow, in section|31 For some aspect ratio, these flow are stable even if they do not fulfill the hypothesis 
of any of the two Arnold's theorems. Applying Arnold's ideas to this case, we first prove their Lyapounov 
stability in section 14.11 In section 14.21 we show the results of direct numerical simulations of the 2D Euler 
equations (nonlinear), which both illustrate the theoretical results of section [51 and show that the linearized 
dynamics correctly describes the nonlinear one. 

Section [S] discusses briefiy some consequences of these results for the 2D Euler and Navier-Stokes equations 
with stochastic forces and for possible theories relying on a quasi-linear or kinetic approach. It also discusses 
some possible generalizations to models of interest for geophysical flows. 



2 The 2D Euler and linearized Euler equations 

Let us consider the 2D Euler equations 

^+V.Vf^^O, (3) 

where O is the vorticity and V is the velocity. We consider this equation either in an infinite plane, in a 
channel geometry with boundary conditions V.n = on the boundary wall, or on a doubly periodic domain 
(0,27r/(5) (0,27r), where (5 > 1 is the aspect ratio. In some parts of the discussion, for technical reasons the 
boundary conditions will be important. Then only the case of a doubly periodic domain will be explicitly 
treated. However all the results are applicable to the channel and infinite domain geometries with slight 
modifications. 

We study the asymptotic stability of parallel flows vo=U {y) e^- We will thus consider the Euler equations 
([3]) with initial conditions close to this base flow : = ljq + lj and V = v + vq, where ojq {y) = —U' (y) is the 
base flow vorticity and uj and v are the perturbation vorticity and velocity, respectively. 

We also need to consider the linearized 2D Euler equations close to this base flow. It reads 

^ + v.Vwo + Vo.Vw = 0. (4) 



We assume that the base flow U {y) has no unstable mode (a precise definition of modes will be given along 
the discussion). In section 2] we will illustrate some of the results on the particular case of the Kolmogorov flow 
U{y) = cosy (in a doubly periodic domain). 

2.1 The Laplace transform, resolvent operator and Rayleigh equation 

In this section and the following we consider the linearized 2D Euler equations. We give the main definitions 
used later on. 

We decompose the perturbation vorticity in Fourier series for the x variable only. For parallel flows, 
due to the translational invariance, these Fourier modes are independent one from the others for the lin- 
ear dynamics. In the following, we thus study perturbations of the form u}{x,y,t) = uJk (y, i) exp (ifcx) and 
ip{x,y,t) = ipk [y T t) exp (ikx) , where tp is the stream function, with w = Atp and k is the longitudinal wave 
number. In the following, we drop the k subscripts for the perturbation. The relations between ui, v and ip are 
then 

.. = ^-fcV, = ^ndvy^^k^. (5) 

dy ay 

The linearized Euler equations then reads 

^+ikU{y)LJ-ik^U"iy)^0 (6) 

We will study the long time asymptotics of the linearized equation. The more general approach is to use 
the Laplace transform of equation Q . We define the Laplace transform a) of w as 

<^{y^p)=l dtuj{y,t)eyip{-pt) 
Jo 
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The Laplace transform is analytic for any complex number p for sufficiently large real part ^p. The inverse 
Laplace transform is given by 

^{y,t) = 7r- dpu}{y,p)exp{pt) (7) 
2vri Jr 

where the complex integration is performed along a Bromwich contour F in the complex plane of p. In the 
following we use the notation p = —ik (c + ie) where c and e are real numbers ; c and e are homogeneous to 
velocities. We assume fc > 0. The Laplace transform uj is thus analytic for sufficiently large e. 
The Laplace transform of equation © reads 

([/(y)-c-*e)^-f/"(y)^=^^^, (8) 

where (j) = tp \s the Laplace transform of and oj (y, 0) is the initial value for the vorticity field. We have 

w = — - k <p 

ay 

From a mathematical point of view, we have to solve the equation for (/) 



with the boundary condition for (p (here doubly periodic). The solution of this boundary value problem 
[aj(., 0)] (y, c + it) depends functionally on the initial value of the vorticity lo (y, 0) (0 is the resolvent operator 
for the stream function). This resolvent operator encodes all the information about the temporal evolution of 
the stream function and the vorticity field. 

The homogeneous part of equation © is the celebrated Rayleigh equation. It is also the equation for modes 
{ip — (p (y) exp [ik (x — (c + ie) t)]} of the linearized Euler equation It reads 

with the flow boundary conditions. For neutral modes (e = 0), this is a non-classical boundary value problem, 
because of the possible singularities associated to the vanishing of J7 — c. Any yc such that U (ye) = c is called 
a critical point for the velocity c. For any c, the free motion on the streamline y = yc is called the critical layer 
and has exactly the frequency kc. 

Any 2/0 such that U'{yo) — (no shear, for instance for velocity extrema) is called a stationary point of the 
jet profile, corresponding to a stationary streamline. We then call c„ — U (yo) a stationary velocity. If yo is a 
local extrema of U, we note that when y — >■ yo (or equivalently c — > co = C/ (yo)), two critical layers, one on 
each side of the velocity extrema, merge on a single one. 

The range of the profile U is the ensemble of velocities c such that miuj, U (y) < c < maxy U (y). 



In the following we assume that the base fiow is spectrally stable, i.e., no unstable mode exist, which means 
that no solution to (fTU| exist for any c + ie with strictly positive e > 0. Then, as shown in section [^751 equation 
© has a unique solution for any c + ie with strictly positive e > 0. We also assume that no neutral mode 
exist, which means that no solutions to (jlOl) is found in the limit e — )■ 0^(a more precise definition will be given 
bellow in terms of the dispersion relation) . This no mode assumption for equation ^ may seem strange, but it 
is indeed a generic situation when such fiows are stables. It is indeed a classical result that shear flow without 
inflection points, or vortex with strictly decreasing vorticity profile are stable and have no neutral mode |23II12| . 
For instance, in section 01 we prove that this hypothesis is also verified for the Kolmogorov flow U [y) = cos (y) 
as soon as the aspect ratio 6 is larger than 1. Actually the only examples we are aware of, of stable flows for the 
2D Euler dynamics, with neutral modes that remain when we add small perturbation to the flow, are cases with 
localized vorticity profile [57] • Usually modes appear, when a parameter is changed, at the edge of an instability 
; four unstable (or two degenerate) unstable modes eigenvalue then emerge from the continuous spectrum. 

The case with neutral modes could be treated following the same lines as what will be discussed bellow ; one 
would then have to separate the contributions by the neutral modes from the contributions by the continuous 
spectrum. The following discussion analyzes the contributions by the continuous spectrum only. 
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2.2 Large time asymptotic for the linearized Euler equation 

In this section, we predict the large time asymptotic of the hnearized Euler equation, using Laplace transform 
tools. We prove results ([T]) and ^ and the mechanism of vorticity depletion at stationary streamlines. The 
heart of the proof relies on the study of the effect of critical layers, on the inhomogeneous Rayleigh equation 
in the limit e goes to zero. This rather technical part is performed in section (|2.3p . 

The results of section (|2.3p are that the resolvent streamfunction (j), solution of has a finite limit for 
small positive e : 

(j){y,c + ie) (j)+{y,c), (11) 

even if singularities exist due to the critical layers. We prove that for any y, (f>^ {y, c) is twice differentiable with 
respect to c, except for velocities c that are in the interior of the range of U . In this last case, for velocities c 
that are not stationary, (/)_|_ is twice differentiable, except for c = U (y). For c — U (y), is continuous but not 
differentiable there, and has a logarithmic singularity: for fixed y, 0+ (y, c) — Aipc {c — U (y)) log {c~U (y)) + 
R (y, c), where R (y, .) is an analytic function of c. When c = cq is a stationary velocity, (y, cq) is differentiable 
with respect to c. 

We think that all of the steps of these proofs could be easily made rigorous from a mathematical point of 
view, by making explicit the required hypothesis. An exception is for the limit of in the case of critical layers 
for stationary points. We actually prove in the following that a limit solution exist for e = 0, but we do not 
prove the convergence to it when e — > 0. In order to deal with this small gap in the proof, we will show, by 
numerically computing (/)(y, e), that this convergence actually takes place. 

These results (the limit and its properties) are the difficult aspects of the discussion, from a mathematical 
point of view. Their technical proof can be skipped at a first reading, the next sections can be read independently 
by assuming these results. The discussion then follows by performing the inverse Laplace transform and proving 
results (H]) and ^ in sections [^3^ and [^331 



2.3 Limit for e — )• 0+ of the resolvent operator 
2.3.1 The dispersion relation 

The equation defining the resolvent operator for the stream function (jH]) is of the type 

dy 



,2 . liy)^ = /(y) (12) 



with q = —k'^ — U" / {U — c — ie) and / = w(y, 0)/[ik{U — c — ie)]. This is a boundary value problem. In order to 
be precise, we treat the case of a doubly periodic domain with the period 2tt, although that is easily generalized 
to the case of a flow in a channel y G (a, b) with the boundary conditions (j){a) = (t){b) = 0. 
For e ^ 0, the differential equation is not singular. We consider the homogeneous equation 

^ + g(y)0-O. (13) 

We consider two independent solutions to p3)) : <j>i is defined by (pi (0) = 1 and (0) = 0, and (j)2 is defined 
by (p2 (0) — and (/)2 — 1 (here and bellow, primes are derivatives with respect to y). The classical variation of 
the parameter computation then insures that a particular solution to (jl2p is 

</'p iy) = -( r ^2f) 4>i iy) +( r 4>if) ^2 (y) , 



and a general solution is 



4>f = 4>p + (14) 



where a and h are unknown constants. The necessary and sufficient conditions for to be periodic are that 
0(0) = (j) (27r) and 0'(O) = 0' {2tt). These conditions read 
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This system has a single solution if and only if the determinant of M is nonzero, which gives the dispersion 
relation 

D{c + It) = [01 (27r) - 1] [0^ (27r) - 1] - <i,\ (27r) 02 (2^) = 0. (16) 

The existence of modes (nontrivial solutions to ([T^ ) is then equivalent to the zero values of the dispersion 
relations. When no mode exist, D is nonzero and (jl2p has thus a unique periodic solution 0j (jl4p . with a and 
h the unique solution to ([T5|). 

Turning back to the inhomogeneous Rayleigh equation, the preceding discussion applies as soon as e 7^ 0. We 
assume that no unstable mode exists, then Z) (c + it) is nonzero. Then the inhomogeneous Rayleigh equation 
has a unique solution for any c + it, for nonzero t. 

The limit e — >■ of 0/ (c + ie) is nontrivial due to the existence of critical layers j/c (c), for which the Rayleigh 
equation becomes singular. We study this limit in the following sections. 



2.3.2 Limit e -H- 0"*" for isolated critical layers 



We consider fixed values of c which are on the range of [/: minj, {U (j/)} < c < maxj, {U (y)}. In such a case, for 
any value of c, there exist one or several points yi such that U{yi) = c. The inhomogeneous Rayleigh equation 
is then singular at such critical layers. We discuss in this section the case U'{yi) ^ (isolated critical layers). 
The case U'{yi) = will be treated in the next section. 

In order to properly study the limit e — >■ 0, we first build a solution to the homogeneous equation pOp . which 
is regular at one of the critical layers y = yi. We define 0^ {y, c) as the solution to ([TU]) with 0^ {yi,c) — and 
4'r iUhc) = 1. From ^TU\\ . we have 0,. (j/;,c) = U"{yi)/U'{yc)- We then have the expansion 



(l^r (y,c) ^ (y-yi (c)) 



u"iyi (c)) 

2U' {yi (c)) 



(y-yi ic)) + o{y-yi) 



(17) 



It can be shown that the solution cjjr {y,c) is an analytic function of y in the vicinity of yi, if we suppose that 
U (y) is analytic in a vicinity of yi. Moreover, from the definition U {yi (c)) = c, because U' (yi) ^ 0, then yi (c) 
is analytic in a vicinity of c and dyi/dc = 1/U' {yi). The solution 0^ {y,c') has then an analytic continuation 
for complex c' in the vicinity of c. 



A classical result of the theory of differential equation of second order is that, if we already know a solution 
(j)r, all other solutions are expressed in terms of 0^ by quadratures. The recipe for this is to look for solutions 
under the form (f) = u(f)r, look for the equation verified by u and integrate it. We apply this recipe to the 
inhomogeneous Rayleigh equation Then any solution to ([H]) is expressed as 

0(y) = dcbr (y) + 0, (y) with (18) 

ft -s r , (2/2, 0)0^(2/2) 

f(y) = / dy2 — , 

Jyo ik{U(y2) - c - It) 

and where d and e are constants. 

We study the behavior of the previous expression close to yi. We first note that / is analytic close to yi. 
Then using the expansion (jl7p . we conclude that 

<p{y,c + it) ^ d(j>r{y,c + it) + g(f>riy,c + it) log(y - yiic + it)) + e(f>g {y,c + it) + (j)h {y,c + it) (19) 

where (jjg and 0^ are analytic functions of y close to yi, and where g is a constant that depends on f{yi), f (yi) 
and e. 

The interpretation of ([T^ depends on which determination of the logarithm we use. Using dyi/dc = 1/U' [yi) 
(discussed above), we have yi (c + it) = yi (c) + it/U' (yi) + o (e). We choose a determination of the logarithm 
such that log {y ~ yi{c) — it/U' [yi)) remains analytic for positive e. Then the study of the limit e — ^ of 
equation (I14p is easily done, we denote this limit (y, c + zO). Using that 0^ and (f'g depends analytically on c, 
we obtain 

0(y, c + iO) = d0r(y, c) + gipriy, c) log |y - y;| + e^g (y, c) + (j)h (y, c) for y > yi and (20) 
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<l){y,c + iO) = {d-iTrsg-n{U'{yi))g)(l)r{y,c)+g(f)r{y,c)\og\y-yi\+e(f>g{y,c) + 4>h{y,c) for y < yi (21) 



where sgn ([/'(?/;)) is the sign of U'{yi). 

From this, we conclude that, for given d and e, the solution to the inhomogeneous Rayleigh equation (|9]) 
converges, for e — > 0, towards a function 0(y,c + iO), which is an analytical function of y, except for y — yi 
where it has a logarithmic singularity. It is continuous at y — yi. This result is valid close to any single critical 
layer yi. If two or several critical layer yu exist in the interval y € (0,27r), then the result is easily extended, 
and holds with a singularity at each critical layer. 

Is this result also true for the solution (j)u:{y, c + ie) of the inhomogeneous Rayleigh equation with boundary 
conditions ? In order to answer this, we now turn again to the construction of section [2.3.11 The result of the 
previous paragraph is applied alternatively to 4>i{y, c + ie), 02 {y, c + ie) and to (j)p{y, c + it). We thus conclude 
that all these three functions have well-defined limits for e — ?> 0+, and that these limits are continuous functions 
of y, which have logarithmic singularities in their derivative for each critical layer. We can then extent the 
definition of the dispersion relation to e — >■ 0+ with _D+ (c) = limj^o+ D {c + it). D+ (c) verifies for which 
we have proved that all terms have a finite limit when e — >■ 0"*". Then we conclude that the two parameters a 
and 6 in equation ()14p have finite limits when e — > 0^. These limit values are given by equation (|15p. where 
each term has a finite limit. 



We thus conclude that the solution to the Rayleigh equation with boundary conditions (j)i^[y,c + ie) has 
a finite limit (j/,c) for e — >■ 0+. Moreover, (y,c) is a continuous function of y that has a logarithmic 
singularity at each critical layer, giving a finite jump for the first derivative. 

Let us denote this jump. From the previous analysis we know that 

0+ {y,c)^a + h{y- yi) + A0+(y - yi) log \y - yi\ + O {y - yif (22) 
Using this expansion, a direct analysis of the dominant term in ^ leads to 

uj{yuQ) + ikU" {yi)<P+ {yi,c) 



tkU' (yi) 



The jump in the derivative thus depends on the value of which is a non local quantity (0+ depends on the 
whole profile U). 

Because (jjr and yi are analytic functions of c, the construction of can be extended analytically when c is 
varied. Then from ((22|) . using yi (c') ~ yi (c) + (c' — c)/U'{yi), one sees that, for fixed y: 

0+ (y, c) = ([/ (y) - c) log {U (y) - c) + 0, (y, c) , (23) 

where (pa {y, c) is analytic close to c = U (y) and with 

^ ^ = Jh^' ■ ^ ^ 



We illustrate the preceding results with numerical solutions of the inhomogeneous Rayleigh equations on a 
doubly periodic domain, and the base flow U (y) = cosy (the Kolmogorov flow). 

We follow the algorithm described in section [2.3.11 that is, computing (f>i, (f>2, and (f>p, and then by using 
them, computing the solution to the inhomogeneous Rayleigh equation for c' = c + ie for small but nonzero 
value of e. In order to numerically compute the solution to the differential equations (for 01, 02, and 0p), we 
use an adaptive method to deal with the singularity close to the critical layers. An extreme precision is required 
in order to obtain satisfactory results. 

In order to test the quality of the numerical simulations, we compute the Wronskian W — (y) 02 (y) — 
02 (y) 01 (y)- From the general theory of differential equations of second order, we know that W does not depend 
on y. Here, from the values of 0i and 02 at y = 0, given by their definition, we have W = 1. We test the 
accuracy of this in all our numerical simulations. For instance in the case of simple critical layers, using the 
Matlab function ode^S, and fixing the relative error and the absolute error parameters of this function to 10""'^^, 
we obtain solutions for which errors in W are typically smaller than 10^^ for e — 10^''. 
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Figure 2: The real upper set of curves) and the imaginary part (G, lower set of curves) for the solution 
4>u) {y,c + i€) to the inhomogeneous Rayleigh equation (0), in the case uj{y,0) — I, c = y/2/2 {yi = 7r/4 and 
yi = 77r/4), k = 1.5. The different curves show the results for e = 10"^(blue (3?) and Hght blue (3)), e = 10"^ 
(green (3?) and magenta (9)) and e = 10"'' (red (3?) and yellow (5)). The curves for e = 10~^ and e = 10"'* 
are nearly indistinguishable, showing good convergence. 
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Figure 3: Same as the previous figure, but for the imaginary part of the derivative 0^ (y, c + ie) of the solution 
to the inhomogeneous Rayleigh equation ^ (blue for e = 10^^, green for e — 10^^) and red for e ~ 10^'^. The 
green and red curves are nearly indistinguishable. 



Figure [2] shows the real and the imaginary part for the solution <j>^ {y, c + ie) to the inhomogeneous Rayleigh 
equation in the case u}{y,0) — 1, c — {yi — n/A and yi — 77r/4), k = 1.5 for the values e — 10~^, 

e = 10^'^ and e = 10~^. This illustrates the convergence of the solutions (pi^ {y, c + ie) to a continuous function 
4'ui {y,c + iO). The visible kinks close to the critical layers suggest the discontinuity of the derivative. This is 
actually verified and illustrated in figure [3J that shows the derivative with respect to y, ipL iu^ c + ie). 



2.3.3 Limit e — 0+ for critical layer of stationary streamlines 

We now consider the case of a critical layer that correspond to a stationary streamline {yi = yo, where U (j/o) = cq, 
t/'(yo) = Oand U" (2/0)^0). 

In order to properly study the limit e ^ in this case, we first build a solution to the homogeneous equation 
(|10p . which is regular at the critical layer y — yo. We define (j)r {y, cq) as the solution to (|10l) with (j>r {yo, Co) = 0, 
(j)^ (yojCo) = and (j)^ (2/0 iCq) = 1. Such a solution can be shown to exist by a series expansion in powers of 
(y — yo). It can be shown that 4>r (y, cq) is an analytic function of y in the vicinity of ?/o, if we suppose that 
U (y) is analytic in the vicinity of yo. 

However, by contrast to the case of isolated critical layers analyzed in section (j2.3.2p . the solution (pr {y, cq) 
is not analytic in the vicinity of cq . Then the approach of the previous section can not be generalized. 

Let us first prove that it exists a solution to the inhomogeneous Rayleigh equation, for e = 0, which is 
continuously differentiable at yo. We start from expression / is analytic in yo. Let us choose e = — /(j/o) 

and 6 = 0. Then we obtain the particular solution 

^..oiy)^My) r^-^^^r^- (25) 

Jya 

Noting that the expansion of ipr begins at order 2 in (y — yo), we easily prove that 

<l>ufi{y) = (y) + 9<f>r (y) log \y-yo\, (26) 
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where 5 is a constant and (pi is an analytic function of y. We note that (pufliu) is continuously differentiable at 
the critical layer t/Q. This solution is defined locally, in an interval where pr has no other zero than yo. However, 
it can be extended to the whole interval y E (0, 27r), because equation © is not singular in other points than 
yo- 

By contrast to the situation obtained for isolated critical points, we can not make an analytical continuation 
of the solution (pst for complex cq + ie. We thus follow another route. 

We note that we can add bpi to 4>u>.a{y)^ where two different values for b can be chosen for y < yo and for 
y > yo- The function 

'/'+(y,co) = b-(pr{y) + PuAv) for V <yo 
'l}+iy,co) = b^ (j)riy) + (jjojAy) y > yo 

is actually a solution to the inhomogeneous Rayleigh equation (|9]) for any y ^ yo which is continuously differ- 
entiable in yi. It is thus a solution to 

The values of 6+ and b^ can be chosen in order to satisfy the boundary conditions. For instance, for 
27r-periodic solutions, the boundary conditions are equivalent to 

The determinant of N then plays the role of a dispersion relation for neutral mode associated to the stationary 
streamlines. It reads 

Ds - (0) (27r) + Pr (27r) (0) (29) 

When no such mode exists, equation (pS)) is solved, and we obtain a solution to the inhomogeneous Rayleigh 
equation that verifies the boundary conditions. 

We have thus constructed a solution to the inhomogeneous Rayleigh equation for real cq ~ U (yo), where yo 
is a stationary point of U such that U" {yo) ^ 0. 

We illustrate the preceding results with numerical solutions of the inhomogeneous Rayleigh equations on a 
doubly periodic domain, for the Kolmogorov base flow U {y) — cosy. 

The numerical computation follows the same rules are the one described in section 12.3.21 We note that 
using the Matlab function odeJ^S, and fixing the relative error and the absolute error parameters of this function 
to 10"""^^, we obtain solutions for which errors in W are typically smaller than 10^^ for e = 10"'^. It is thus 
much harder to obtain good quality numerical simulation in that case, than in the case of isolated critical layers 
discussed in section [2.3.21 

Figured] shows the real and the imaginary part for the solution <j)i^ (y, co + ie) to the inhomogeneous Rayleigh 
equation in the case 0) — 1; with a critical layer corresponding to a stationary point (cq = 1, j/o = 0) 
and with k — 1.5 ; for the values e = 10^^, e — 5.10"^ and e = 10"'^. This illustrates the convergence of the 
solutions (pi^ iy,co + ie) to a continuous function {y,CQ). It turns out that the real part converges to zero. 
The same results are also presented for the derivative with respect to y, p'^ (y, cq + ie), on figure [5j 

We now turn to the derivation of a property of such solutions, that will be very important in the discussion 
of the asymptotic behavior of the linearized 2D Euler equations. We have shown that {y, cq) is continuously 
differentiable at yi, and has a second order logarithmic singularity at yo (see ((26)) and ([27l) l. Then a direct 
inspection of the leading singularity in equation of order (y — yc)~^, leads to the conclusion that 

tkU" (yo) <P+ [yo, U (yo)) + c^(yo, 0) = (30) 
2.3.4 The asymptotic vorticity field 

Using the results of the previous section, we prove in this section that the vorticity field converges, for large 
time, towards a field oscillating at a multiple of the streamline frequency. More precisely, we prove that 

uj{y,t)^ Lo^iy)e^pi-ikUiy)t) + o(^) (31) 

In particular, for any stationary point yo, uJoo (yo) = 0. This essential property means that the vorticity cancels 
rapidly at any stationary streamline. This is the mechanism of vorticity depletion at the stationary streamlines, 
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Figure 4: The real (3?) and the imaginary (9) part for the solution (p^ (y, cq + ie) to the inhomogeneous Rayleigh 
equation ([9]), in the case cij(2/, 0) — 1; with a critical layer corresponding to a stationary point (cq = I, yi — 0) 
and with k = 1.5. The different curves show the results for e = 10~^ (blue (3?) and light blue (3)), e = 5.10"'^ 
(green (3?) and magenta (5)) and e = 10~^ (red (3?) and yellow (5)). The curves for e = 5.10"'^ and e = 10"'^ 
are indistinguishable, showing good convergence. The real part is the ensemble of curves that converge to zero. 




Figure 5: Same as the previous figure, but for the imaginary part of the derivative (j)'^ (y, cq + ie) of the solution 
to the inhomogeneous Rayleigh equation All three curves are indistinguishable. 
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discussed in the introduction. 



Using (fTTj) and ([8]) we have 



6^0+ ifc (f7(y) - c - iO+) ^ ^ 



Thanks to the analysis of the properties of (f> in the previous section, we know that all its singularities are 
integrable. We thus see that for any c, w(.,c) has non integrable singularities at each critical layer yi. For fixed 
y, Cj{y, .) has a single singularity for the velocity c = U{y). Using (|32p . we write the inverse Laplace transform 
([7]) on the Bromwich contour defined hy p = —ik{c + ie) with e > and — oo < c < +cxd. 

We first estimate the contribution of the pole 

1 /•+°° exp Hfc(c + »e)t) j^^^„ C/(y))+a.(0, y)] = [zfcC/" (y) 0+ (y, [/(y)) + w (y, 0)] exp {-zkU{y)t) , 

2m J _^ U[y)-c~ie 

using the standard deformation of the contour of the complex integral, and the residue theorem. The remainder 
contribution to the vorticity (|33p is then 

^"°°dc^(^iii^|±(^i«exp(-..ci). 



2^ 7-00 U{y)-c 

This integral is an oscillating integral. For large times, it thus gives a contribution of order O (l/i'^) where 7 
depends on the order differentiability of 0+ (y, c) as a function of c. The result (j3ip is thus proved, and we have 

Woo (y) = ^kU" (y) (y, U{y)) + w (y, 0) . (34) 

We remark that for any point yi where U"{yi) vanishes, a;oo(yi) = w(yi,0). This could have been antici- 
pated as for such points yi, from ([5]), we trivially have lo (yi, t) — w(yi, 0) exp {—ikU (yi) t) for any time t. 

Using and we deduce that 

Woo (yo) = 0. 

This result means that the vorticity tends to zero for large time for any stationary streamlines yp. This vorticity 
depletion at the stationary streamline is, from a mathematical point of view, a nontrivial consequence of the 
Laplace transform analysis, and of the regularization of the resolvent operator at stationary velocities. As will be 
illustrated in section using numerical simulation, it is a striking dynamical effect leading to the disappearance 
of any filament in the area close to the critical layer of a stationary point of the profile U . This has a large 
qualitative impact on the flow structure and evolution. 

This effect comes from the term VyU" (y) in the linearized Euler equation ([H]) ; it is thus a consequence of 
the effect of the transverse velocity on the background vorticity. Because Vy is a non local quantity, depending 
on the evolution of the vorticity field everywhere in the domain, this effect is a non local, non trivial one that 
we are not able to explain easily heuristically. 

Besides these theoretical results, the Laplace tools are very interesting as they allow the computation of 
asymptotic behavior of the flow without relying on a complex direct numerical computation. Moreover, whereas 
in asymptotic approaches like the Lundgren's one, where the asymptotic profile Wqo (y) is not determined, here 
we can compute it from (j34p . 

Using this last procedure and the numerical computations of the resolvent (pi^ , described in sections 12.3.21 
and 12.3.31 we compute the asymptotic vorticity profiles. They are represented in figures [TUl [TT] and [1^] ; and 
discussed in more details in section [ 
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2.3.5 The asymptotic velocity field 

In this section we study the asymptotic behavior of the velocity field. We prove that the velocity field decays 
algebraically for large times: 

, ^oo jy) exp {-ikU{y)t) 



^oo jy) exp{-ikU {y)t) _ 
tk{U'{y)f 

where Woo is the asymptotic vorticity profile pil34p . 

We first explain this result starting from the asymptotic vorticity derived in the previous section pip , and 
using large time asymptotic behavior of oscillating integrals. This argument is heuristically very interesting. 
However it will be valid only when the contributions of stationary points are negligible, and when the 
convergence of the vorticity towards the asymptotic vorticity is sufficiently rapid. Indeed, the convergence 
towards the Lundgren's ansatz has to be uniformly more rapid than the derived algebraic laws. This last point 
is actually true only for strictly monotonic velocity profiles U as will be discussed below. 

In order to give a proof of the results ((35)) and (l35t valid also for profile U with stationary streamlines, we 
give a more general argument based on Laplace transform at the end of this section. 

Oscillating integrals We begin with the expression of the velocity from the vorticity field using a Green 
function formalism. We have 



v(2/,t) = j dy'Gk{y,y')u:{y',t), (37) 

where is defined from ([SJl: 

Gk{y,y') = [-^,ikHkyy,y') with ^ ^ k^H, ^ S (y - y') , 
with the periodic boundary conditions on y. Using the asymptotic result on the vorticity field, we thus have 

v{y,t)^^ j dy'Gk{y,y')LOoo{y')e^Y>{-ikU{y)t). (38) 



We consider the asymptotic behavior, for large times i, of the oscillating integral ([55)1 . Since Kelvin, very 
classical results do exist for the asymptotic behavior of such integrals, the most famous result being the stationary 
phase approximation. Such results are discussed in appendix A. 

An essential point, which makes this case different from the more classical ones, is that the Green function 
Gfc (y, y') is not smooth everywhere: it is smooth except for the singularity when y = y' . We prove in appendix A 
that ifU (y) has no stationary points (for all y, U' (y) ^ 0), then results (|35ll36p are valid, the main contribution 
being related to the singularities of the Green function. 

If the velocity field U (y) has M stationary points t/m {U' [ym) ~ 0), then the contributions of the stationary 
points, generically of order l/Vt, usually dominate the contributions of the singularity of the Green function 
Gfc, in integrals like ()38p . If {ym) 7^ 0, the classical stationary phase approximation (see appendix A) would 
lead to 



: yrn)^oo {yn 



27r / iemTT\ exp(— i/ct/ {ym)t) 



m=l..M 



fct/"(2/™)i ''''p V 4 ; vt 



where em is the sign of —kU" (ym)- 

However, a remarkable fact is that for any stationary point j/o (such that U' (ya) = 0), due to the vorticity 
depletion mechanism discussed in section 12.3.41 and proved in section 12.3.31 Woo (yo) = 0. Then the leading 
order contribution from the stationary phase approximation vanishes. The analysis could proceed in order to 
determine the next leading order term in the expansion, from ()38p . expected to be of order l/t^^"^. However, 
such a detailed analysis is useless, because the convergence of a;(., t) towards the asymptotic vorticity profile Woo 
is too slow, in the vicinity of a stationary streamline. Actually the error due to the slow convergence towards 
the Lundgren's profile gives contributions which are also of order This will be illustrated using direct 

numerical simulations in section |3|(see figure [Hj page [M) figure [T5] page HUl and the related text). 
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Laplace tools In order to give a precise argument for the results (I35ll36p . we use Laplace tools. We first 
note that Vy = ikip, and study the asymptotics for the stream function ip. Starting from the inverse Laplace 
transform of ip, we have 

V'(y,i) = ^/ dc(j)+{y,c)exp{-ikct), (39) 



where 0+ (y, .) is the limit of (/) {y,c + ie) for e — >■ 0. (j39p is an oscillating integral. We use that for any y 
for which U' [y) ^ 0, 0+ (y, c) is twice differentiable except a± c — U (y) where it has a logarithmic singularity 
A(j)c {c — U (y)) log {c — U (t/)) (see equation ([121) )• Then the large time asymptotics of -0 is due to this singularity. 
In order to evaluate it, we part integrate twice (j39p and evaluate the contribution of the singularity with the 
residue theorem. Then the leading order contribution is obtained as, 

{tkU'iy)f t2 ' ^ ^ 

where we have used and in order to express A0c- 

We note that (j36p follows immediately from (PO)) and from the relation Vy = ikip. The asymptotic result 
(j35p for the transverse velocity can be derived by following similar arguments as the one just described for 
0. 

The above argument uses the explicit prediction (j23p for the singularity of 0+ (y, c) as a function of c. The 
expressions and ([^^ are valid only when y is not a stationary streamline ([/' (y) ^ 0). For stationary stream- 
lines j/m we have no theoretical predictions. Direct numerical computation in section l¥^ will lead us to conjecture 
that for such special points i/j (j/m, ^ Ciex.p {—ikU{ym)t) /t^^'^ , Vx {ym,t) ~ C2e:xjp {—ikU{ym)t) /t^^"^ 

i— >C30 i— >-C30 

and ti„ (ym, t) ^ ikCi exp l—ikU {ym)t) /t^/"^ . We note that this exponent 3/2 is not related to a contribution 
from the stationary phase approximation. 

We thus conclude that the results (|35p and (|36p are valid for monotonic profiles, but also for 
base flow with stationary streamlines. This is in marked contrast to what was thought in many 
previous publications based on the asymptotic expansions and the stationary phase arguments. 

This is mainly due to the vorticity depletion mechanism at the stationary streamlines discussed in the previous 
section. We also stress that, using Laplace tools, the asymptotic profile Voo [y) can be numerically computed 
easily, without relying on direct numerical computations of the Euler equations. 

We have theoretical predictions for the power law in the asymptotic behavior of the perturbation velocity, 
for all points of the domain except along the stationary streamlines. 



3 Asymptotic stability of parallel flows for the Euler equations 

In the previous section, we have obtained results for the asymptotic behavior of the linearized Euler equations, 
with initial conditions close to some parallel flows vq (r) = U [y) e.^ . We now address the evolution of the same 
initial conditions by the nonlinear Euler equation ([3]). The aim of this section is to explain why the linearized 
dynamics will be a good approximation for the dynamics for any time i, and to explain why the flow velocity is 
asymptotically stable (in kinetic energy norm), for small initial perturbation of the vorticity (in the enstrophy 
norm). 

We consider the initial vorticity (a;, y, 0) = —U' (y) -few (x, y, 0), where e is small. We suppose, without loss 
of generality that J dxoj = Q. The perturbation uj can be decomposed in Fourier modes along the x direction 

w {x, y,t)^^ ujk (y, t) exp [ikx) . 

k 

From the Euler equations the equation for the evolution of uj^ reads 

^^ikU{y)iOk~-iki^kU" (y) = ~eNL with NL = Y, {"*^^ iy,t)^k~i {y,t) + ^ [il^Ji {y,t)ujk-i (y,t)]| . 

(41) 

The left hand side is the linearized Euler equation, whereas the right hand side are the nonlinear corrections. 
We want to prove that, for sufficiently small e, neglecting the nonlinear terms is self-consistent. 
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For this we have to prove that the nonhnear terms remain uniformly neghgible for large time. We then use 
the asymptotic results for the linearized equation derived in section We thus have, for any k 

A,L{y,t) ' ,,2 75 and Wfc,L (y, i) - ujk.L.o^ (y) exp {-ikU{y)t) , (42) 

{ikU' [y)) i 

where the subscript L refers to the evolution according to the linearized dynamics. We call a quasilinear 
approximation to the right hand side of equation (j4ip , the approximation where -0^ and Wfc would be evaluated 
according to their linearized evolution close to the base flow U (y) . From (|^^ , one would expect at first sight 
that this quasilinear approximation of the nonlinear term NLql, would give contributions of order 0{l/t). 
The detailed computation, easily performed from (|32]), actually shows that the contributions of order O (1/i) 
identically vanish for large times. Then 

eNLk^Q, =^ O (±) 

This is an important remark, as it proves that within a quasilinear approximation, the contribution of the 
nonlinear terms NLql remains uniformly bounded, and more importantly it is integrable with respect to time. 

Then it is natural to conjecture that the contribution of the nonlinear terms remains always negligible. More 
precisely, we naturally conjecture, that within the fully nonlinear equation, for sufficiently small e: 

ipk[y,t) n 75 and cjfc y,t - Wfc,oo y exp -ifcC/ y)t 

{ikU' {y)) i 

with 

wfc^oo {y) = wfc,L,oo (y) + O (e) 

A similar reasoning in order to evaluate the nonlinear evolution for the profile U (y) would lead to the 
conclusion that for large times 

no (y, i) , ^ ^U'^ (y) with U'^ (y) = U{y) + SU (y) , 

where SU = O (e^) . 

This means that the parallel flow will quickly stabilize again towards another parallel flow which is close to 
the initial one. 



A natural question would be to compute the modified profile. The preceding analysis leads to the quasi-linear 
expression 

/>oo 

SU (y) = -e^ / dt NLo.QL [t) + o (e^) . (43) 
Jo 

This expression involves integrals over times of the linearized Euler equation. It is not amenable to a simple 
expression, but could be easily be computed numerically from Laplace tools. 

This result has to be contrasted with the results usually obtained using a quasi-linear approach, for instance 
in the kinetic theory of particle dynamics (point vortex models, plasma physics, astrophysics). Usually the 
integrals occurring in (|43p diverge. Then one invokes a time scale separation, and the divergence of the integral 
is regularized using a multiple scale analysis. Here, by contrast, the integral converges. This means that there 
is a single time scale over which all quantities reach their asymptotic value (a typical time scale here is 1/s 
where s is the typical shear). The nonlinear evolution is thus very brief and leads to very small changes in the 
initial profile U . 

A theory for the relaxation towards equilibrium for the 2D Euler equations has been proposed based on 
a quasi-linear theory coupled to some Markovianization hypothesis, by analogy with the kinetic theory of 
point vortices [501 ■ Such an approach based on analogies are natural guess that rely on theoretical hypothesis 
(quasi-linear hypothesis and Markovian hypothesis), that would benefit from either theoretical justifications or 
numerical verifications. The results of this paper for the linearized dynamics and the discussion above, show 
that the approximation of the non-linear dynamics by the linearized dynamics remains uniformly self consistent, 
and a quasi-linear approach self-consistent. However as discussed above, the relaxation is then extremely rapid 
and all quantities relax rapidly. Any further assumptions like Markovianization seems then irrelevant. Simpler 
approach, like the one of [H], similar to the discussion of this section seems thus more relevant. 

We thus conclude that a direct asymptotic expansion based on a quasi-linear approach, similar to the one 
in |41| or may be with more subtle treatments of the lower order contributions, is most probably the relevant 
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approach for a theory of the relaxation of the 2D-Euler equations. The exact results on the linearized dynamics 
of the previous section and the discussion above put such an asymptotic expansion on a more rigorous basis, by 
proving that the leading order contribution remains self-consistently bounded for all times, and explaining why 
an asymptotic expansion possibly converges. In the Lundgren approach, the asymptotic behavior is described 
by an asymptotic expansion for large times, whose leading order term is not determined. It actually depends 
on the initial condition and can not be predicted only with a large time asymptotic expansion. The further 
interest of the Laplace method developed in this paper is to give precise predictions for the asymptotic profile 
that can be easily computed numerically from the Laplace tool. The numerical computations of next section 
will confirm the statements on this last paragraph, by showing excellent agreements between direct numerical 
simulation of the Euler equations on one hand and the prediction of such a simple quasi-linear approach, based 
on Laplace tools computations. 

We also conclude that, for any profile U verifying the hypothesis of this work (no unstable and no neutral 
modes for the linearized dynamics) with any perturbation corresponding to a small vorticity, the assumption 
that the velocity converges for large times towards a new parallel velocity profile close to the initial profile U is 
a self-consistent hypothesis. We see in section 2] that this is confirmed by numerical computations. 

From this discussion, we thus conclude that it is natural to conjecture that for any profile U verifying the 
hypothesis of this work (no unstable and no neutral modes for the linearized dynamics), for any perturbation 
corresponding to a small vorticity, the velocity converges for large times towards a new paralleled velocity profile 
close to the initial profile U. A possible theorem expressing this more precisely would be similar to the one 
recently obtained by Mouhot and Villani [35], for the Landau damping in the very close setup of the Vlasov 
equation. It has however to be noted that no proof of such a theorem for the Euler equations is known yet, 
even in the simplest case of a profile U without stationary points. 

It is thus very natural to conjecture that the ensemble of shear flows without unstable or neutral modes is 
asymptotically stable in the sense given previously (initial perturbation controlled by a vorticity norm, for 
instance the enstrophy and large time perturbation controlled in kinetic energy norm) Q 



4 The Kolmogorov flow 

In this section, we consider the particular case of the 2D Euler equations in a doubly periodic domain 2? = 
[0, 2Tr/5) [0, 2tt), where (5 > 1 is the aspect ratio ; with the Kolmogorov base flow U {y) = cos (y). 

4.1 Stability 

In this section we study the stability of the Kolmogorov flow for 2D Euler equation dynamics. We note that, 
in the linearized 2D Navier-Stokes equations, Mishalkin and Sinai |35] have found that the Kolmogorov flow is 
stable for 6 > 1 and unstable for 5 < 1. Here we prove the stability for (5 > 1 for the dynamics of both the Euler 
and linearized Euler equations. We also show numerically that unstable modes exist for S < 1, while no stable 
mode exists for 5 > 1 for the linearized 2D Euler equations. 

4.1.1 Lyapounov stability 

We consider initial vorticity conditions close to the base flow vorticity wq (t;) — sin(i;). We let this initial 
condition evolve according to the nonlinear Euler equations If the perturbation to the initial flow remains 
small for all times, the flow is said to be Lyapounov stable. We first prove that the base flow ujq {y) is Lyapounov 
stable as soon as (5 > 1. 

*We refer here to the notion of asymptotic stability of an ensemble of steady states of an infinite dimensional Hamiltonian 
equations, see for example the work [52^ where any stable soliton of the KDV equations, slightly perturbed, is proved to converge 
for large times towards another slightly different soliton. Asymptotic stability of ensemble of steady states has also been proved for 
other solutions of infinite dimensional Hamiltonian systems. 

classical argument, presented in a rigorous framework by Caglioti and Maffei 1171 in the context of the Vlasov equation, 
implies that steady states of the Vlasov equation for which Landau damping would occur, would be unstable in a weak norm. 
At the core of the argument lies the time reversal symmetry of the equations. These arguments would be easily generalized to 
the Euler equations. This may seem in contradiction with the notion of asymptotic stability discussed here. However the notion 
of stability discussed by Caglioti and Maffei involves weak topology for both the initial conditions and final state. There is no 
contradiction with our definition of asymptotic stability, as we control here the initial perturbation in a vorticity norm and control 
the convergence in a velocity norm. 
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The classical Energy-Casimir method proposed by Arnold [2l cannot be applied directly here. Indeed, the 
Kolmogorov flow does not verify the hypothesis for any of the two Arnold theorems [2]. However, we can still 
prove the stability in this case, by direct analysis. Let us define the energy- Casimir functional F as 




F being half the enstrophy minus the kinetic energy, it is a conserved quantity for the Euler equations. 

We first prove that the base flow loq (y) = sin (y) (vq = cos{y)ex) is a minimum of F. We consider the 
perturbation vorticity uj = Q — wq, and decompose it into Fourier modes lo = X)fe>i'-^fe6fc with Ae^ = — AfcCfc, 
where the Afc > are arranged in increasing order, and where Cfc are orthonormal (Jj^ekSk' = Skk')- On the 
doubly periodic domain D = [0, 27r/S) [0, 2tt), if (5 > 1, we have Ai = A2 = 1, corresponding for instance to the 
modes cos(?;) and sin (y). Then for any /c > 3, A^ > 1. 

We obtain 

k 

Since Afc > 1, F > 0. Moreover F [ujo] — 0. We thus conclude that uto is a global minimum for F. 

We note that this minimum is degenerate, as all vorticity fields co = a cos (y) + /? sin (y) are also minima. 
Since _F is a conserved quantity, we conclude that 

^Y^hh_l^l^t)^e,, (45) 

fc>3 

where ep = F (0) is the small value of F for the initial perturbation. Then if they are initially small, all Uk for 
A: > 3 remain small for large times, the amplitude being measured according to the norm (|44p . 

Expression does not control the first Fourier modes fii and f22. For this, we use the fact that the 
enstrophy 

r2= [ 

Jv 

is conserved. We suppose that its initial value is r2.o + er where r2,o is the base flow enstrophy and er is the 
perturbation enstrophy. Using the enstrophy conservation we have 

k 

Then, using that ^^/^>3 ^ x^^—i ^ ^ (derived from (|45p . using X]^ ^ A3 for A; ^ 3), we have 

\Vll{t) + nl{t) - r2,o| < max jer, J-^^f^ ■ 

This means that the flow associated to the two first mode is a{t) sin(?/ + where (j) may be arbitrary but 
where the amplitude a is controlled up to an error of order max{er, ^f}- 

We have thus proved that any initial condition close to the initial profile oj = sin (y) remains close to the 
family of profiles sin (y + 0). Then the flow is Lyapounov stable in this sense. 

4.1.2 Linear and spectral stability 

Next, we let the initial conditions close to the base flow a;o(y) = siny evolve according to the linearized 2D 
Euler equations If the perturbation to the initial flow remains small for this dynamics, the flow is said to 
be linearly stable. 

We decompose the perturbation vorticity in Fourier series for the x variable only. For parallel flows, due 
to the translational invariance, such Fourier modes are independent. The modes with no dependence on x are 
easily shown to be neutral. Then the only issue is about the stability of other modes. In order to prove this, 
we simply note that the perturbed Energy-Casimir functional (I44p is conserved not only by the nonlinear Euler 
dynamics but also by the linearized one. Then, because the Energy-Casimir functional is positive, this proves 
that any a;-dependent perturbation remains small if it is initially small. The flow is thus linearly stable as soon 
as 5 > 1. 

If the linear equation has no exponentially growing modes, it is called spectrally stable. Linear stability 
implies spectral stability (the converse may be wrong) . Then because it is linearly stable, we can conclude that 
no unstable modes exist to the linearized 2D Euler equation as soon as (5 > 1. 
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Figure 6: Minimum values for the dispersions relation mine (c, k) as a function of . This plot shows that 
neutral modes exist only for — 1. 

4.1.3 Neutral modes 

We look for the modes of equation (|4]) such that the stream function is of the form V' = (y) exp {ik {x — ct)). 
Then (p satisfies the classical Rayleigh equation (|10l) . 

As discussed in section |4.1.2[ no unstable eigenvalue exists for fc^ > 1. Thus only for real values of c can 
the Rayleigh equation have solutions for k^ > 1. In the following, using numerical simulations, we show that 
no neutral modes exist, except for the marginal case fc = 1. 

When c is in the range of [/: — 1 = miiiy {U (y)} < c < 1 = maxj, {[/ (y)}, U ~ c vanishes at the two 
critical layers defined by U {yi-^ 2) — cos {yi-^ 2) = c. The Rayleigh equation then has logarithmic singularities. 
As discussed in section 12. 1[ when initial value problems are considered, the relevant solutions to the Rayleigh 
equation are the ones that are obtained with c' — c -\- it, c real, and by considering the limit e 0+. We study 
the existence of modes in that limit. 

For this, we numerically compute the dispersion relation _D_|_ (c, k) of the Rayleigh equation, as defined in 
section [2.3.21 Neutral modes correspond to zeros of We use the same numerical tools as the one described 
in the end of section 12.3.21 we use the Matlab function ode4-5, and fix the relative error parameter and the 
absolute error parameter of this function to 10^^'^, then obtain solutions for which errors in the Wronskian 
W are typically smaller than 10^^ for e = 10^'*. We approximate (c, k) by the numerically computing 
D{c + ie,k) with e = 10""^. 

Figure [5] shows D„i (k) ~ miuc (c, fc) (we note that _D+ is unchanged when the sign of k is changed). For 
a given value of fc, some neutral mode exist if and only if Dm (k) vanishes. We conclude from this plot that 
neutral modes exist only for the value k^ = 1 (we have tested values of fc^ up to fc^ = 10). 

For fc^ = 1, we see numerically that a mode exists for c = only. The mode is then found by direct inte- 
gration of equation pop . It is the trivial mode (p = C where C is any constant number {ip = C exp {ikx) and 
ip = C exp (—ikx)). 

The representation of the complex curve (c) when c is varied is called a Nyquist plot (see [48] or [5] in 
the context of fluid dynamics) . It is very useful, as the algebraic number of loops of the Nyquist plot around 
counts the number of unstable modes on the complex half plane c — cii + iX with positive A |48l [5]. 
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Figure 7: Nyquist plots (complex _D+ (c) represented in the complex plane, when c is varied), for values fc^ — 0.99 
(green dashed line) and fc^ — 1.01 (plain blue line). 

Figure [7] shows the Nyquist plot of £)+ (c), for k — 0.99 and k = 1.01 respectively. One clearly sees the 
passing of the curves across the value — when k is changed from k — 0.99 to k = 1.01, corresponding to 
neutral modes for fc = 1. Moreover, we conclude that only one neutral mode exists for this value of k, because 
only one branch of the curve passes through 0. For larger value of c (not shown), the upper part of the plot 
loops to the right on the upper half plane goes down to the lower half plane by crossing the real axis for very 
large values of (f>, before to close on the branch visible in the lower half plane. Counting the algebraic number 
of loops around zero, we thus conclude that for /c > 1, no unstable modes exists, in accordance with the result 
of section (|4.ip ; whereas for /c < 1 only one unstable mode exists. 

From this analysis, we thus conclude that only one neutral mode exists. This modes corresponds to the 
destabilization of the flow and the appearance of an unstable mode, when passing from values k^ > 1 to values 
k^ < 1. It is the trivial mode tp = C exp (is), found for the phase speed c = 0. 

4.2 Direct numerical computation of the 2D-Euler equations 

In this section, we illustrate and complement the above results through the direct numerical simulation of the 
2D (nonlinear) Euler equations, in the doubly periodic domain T) = [Q,2tt/5) [0,27r), for the Kolmogorov base 
flow U{y) — cos y. 

Since the base flow is homogeneous in the x direction, the dynamics of the fields would be decoupled into that 
of the components of the form fk{y^ t) exp (ikx), if the equations were linearized. Then it is natural to consider 
fields of this form also in the (nonlinear) Euler equations, because we are interested in slightly perturbed parallel 
fiows. 

Since, in the Euler equations, instabilities are mainly large scale ones, the flow is more likely affected by large 
scale perturbations. Moreover, in geophysics and experiments, forcing and perturbations tend to be dominantly 
effective on the largest scales of the fiow. It is thus natural to study initial perturbations of the form 

uj{x, y, 0) = e^(y) cos Sx (46) 

where k = 6 is the smallest wavenumber (largest scale). We consider the case e ^ 1 ; we use e = 0.01 throughout 
in the following numerical computations. 
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Figure 8: The time series of perturbation velocity components | vs^x {y = 0, | , for the initial perturbation profile 
A{y) = 1 and the aspect ratio 5 = 1.1, with various resolutions of the system {N x N). The simulations blow 
up at i 70 for N — 256 and at t 150 for N — 512. In the following numerical results, we have used high 
resolution enough to confirm that the asymptotic behavior is robustly observed. 

Similarly, we first examine the dynamical response for the same wave number as the initial perturbation, 
namely. 



The analysis of nonlinear effect will be performed at the end of this section. 

Direct numerical simulations. In the following, for the direct numerical simulations, we use the classical 
pseudo-spectral method algorithm [S^, which is the most precise and robust numerical algorithm currently 
known for the Euler and Navier Stokes equations in doubly periodic domains, and indeed regarded as the 
standard method for precise numerical simulations in these cases. 

A standard way to compute numerical solutions to the Euler equations is to integrate directly the truncated 
Euler equations without viscosity [33J . We follow this procedure. The numerical scheme is thus not stable in the 
long run; on account of the lack of explicit dissipation, the simulation blows up when the small scale structures 
become of the same size as the grid scale [33 . As an illustration, figure 8 presents the computation of the first 
mode of the velocity perturbation, computed with three different numerical simulations using 256^, 512^ and 
1024^ Fourier components, respectively. This illustrates the blowup after some times, occurring later on for 
larger resolutions. Moreover this shows that the short time numerical results are stable when the resolution 
is increased. In the following numerical results, we have always confirmed that the resolution of the system 
is high enough, by systematically testing the stability of the results by changing the resolution, such that the 




(47) 
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asymptotic behavior discussed is robustly observed. This procedure for assessing the numerical results for the 
Euler equations without dissipation is also a standard one (see for instance [33]). 

An independent assessment of the quality of the numerical results is provided by figure [TOl discussed later in 
more details. On this figure, the final vorticity profile is computed by two algorithms: the just described direct 
numerical simulations, and the predictions from the Laplace transform (equation (34)). The results from these 
two completely independent algorithms are indistinguishable on figure 1101 The maximum discrepancy between 
the two profiles on this figure is of order 0.0001. 

Asymptotic vorticity profile for even perturbations. First, we consider an initial perturbation where 
A[y) has the same parity as the base fiow. In particular, we examine A{y) = 1. 

The space-time series of |a;5(j/,t)| is shown in Figure [1] page [51 which we have already seen. Initially, it 
rapidly (almost exponentially) relaxes toward the final profile, |w5oo(2/)|; in particular, it relaxes to zero at y = 
and TT (stationary streamlines), whereas it remains constant at ?; = tt/2 and 3tt/2. The rapid relaxation of the 
modulus |a;5(y, t)| is in agreement with the theoretical prediction pip that the Lundgren ansatz is asymptotically 
valid. 




y /2n 

Figure 9: The profiles of perturbation vorticity, \ujs{y,t)\, at several times, for the initial perturbation profile 
A{y) = 1 and aspect ratio 5 = 1.1. A flat region is observed near the stationary streamline y = 0- As time goes 
on, this flat region becomes narrower and narrower, and instead the region with the profile proportional to y^ 
extends towards the stationary streamlines, leading to a parabolic profile in the large time limit. 

After the rapid relaxation, |w5oo(2/)| converges towards the final profile algebraically. In particular, in the 
vicinity of y = and tt (stationary streamlines), it relaxes as slowly as t~^, leading locally to a flat profile (see 
figure in]). However, as time goes on, this flat region becomes narrower and narrower, and instead the region 
with the profile proportional to y^ extends towards the stationary streamlines. This indicates that, in the large 
time limit, the profile is parabolic in the vicinity of the stationary streamlines. 

This also illustrates that the relaxation towards the stationary profile does not converge in a uniform way; 
the process is very slow close to the stationary streamlines whereas it is extremely rapid away from them. 
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Figure 10: The final profile of the modulus of the perturbation vorticity, |cj5oo(y)|, obtained from the direct 
numerical simulation and the theory (equation ([34])), for the initial perturbation profile A{y) = 1 and the aspect 
ratio 5 = 1.1. The two profiles show a very good agreement. 
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Figure 11: The final profile for the modulus of the perturbation vorticity, |a;5oo(y)|, for the initial perturbation 
profile A{y) = 1, and aspect ratios S = 1.1, 1.2 and 1.5, computed from the prediction of the Laplace transform 
tools (equation ([M]) '). The profile shows a bifurcation from single to double peak shapes, when S is increased. 
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(b) 

Figure 12: The space-time series of \vs,xiy, t)\ i^) ^^'^ \^s,y{y, t)\ (b), for the initial perturbation profile A{y) = 1 
and the aspect ratio 5 = 1.1. Both the components relax toward zero, showing the asymptotic stability of the 
Euler equations. 
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Figure 13: The time series of perturbation velocity components \vs,x{y,t)\ (a) and \vs,y{y,t)\ (b) at three 
locations, y = (vicinity of the stationary streamline) (red), y = tt/A (green), and y = it/2 (blue), for the initial 
perturbation profile A{y) 1 and the aspect ratio S = 1.1. We observe the asymptotic forms \vs.x{y,t)\ ^ t^", 
with a — 1, and \vs^y{y,t)\ ^ t~^, with /3 = 2, in accordance with the theory for the asymptotic behavior of the 
velocity (equations and ([55)) 1 
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The width yj of the flat area decreases as ~ whereas the constant value of the vorticity modulus 

in the flat area decreases as 1/t. When computing the velocity from the vorticity (equation ([57]) ). the overall 
effect of this flat area is thus of order Such a contribution is thus of the same order as what would give 

the leading order term of the asymptotic expansion of oscillating integrals, as discussed in the end of section 

In Figure [TUl we compare the final profiles obtained from the Laplace tools (equation ([M]) ) and the direct 
numerical simulations. As shown, the results show a very good agreement. This agreement support both the 
quality of the direct numerical simulations and the results of the computation of the asymptotic profile from 
the Laplace transform tools. 

The computation of the asymptotic profile from the Laplace method is extremely rapid and easy, compared 
with direct numerical simulations. Using this tool, we study some qualitative properties of the asymptotic 
profile. By increasing the aspect ratio (5, we observe a bifurcation from a single- to a double-peaked asymptotic 
vorticity profile (figure [TT|) . The three asymptotic profiles all show the depletion of the vorticity perturbation 
at the stationary streamlines. 

Asymptotic decay of the velocity perturbation for even perturbations. The space-time series of the 
modulus of the perturbation velocity components \vs^xiy,t)\ and \vs,y{y,t)\ are shown in 1121 The relaxation to 
zero of the velocity perturbation illustrates the asymptotic stability of the velocity for the 2D Euler dynamics. 

We investigate the asymptotic behavior of the velocity perturbation more precisely. Figure [T2] shows the 
time series at several positions. As shown, their asymptotic forms are \vsx(jj,t)\ ^ t^", with a — 1, and 
\vsy{y, t)\ ^ t^^ , with (3 — 2. This is in agreement with the theoretical predictions for the asymptotic behavior 
of the velocity perturbation (see equations 1551 and 1551) ) . 

Odd perturbations. Next, we consider initial perturbations where A{y) has a parity opposite to the base 
flow one. In particular, we examine A(y) — siny. 

The space-time series of \uJs{y, t)\ is shown in Figure [TD It shows an initial rapid relaxation toward the final 
profile, as expected from the theory (equation ([M)) ). Since the parity of the perturbation is conserved for all 
time, the vorticity profile remains odd. Then |w5(y,i)| is zero for y = and y = tt (stationary streamlines), as 
expected. 

The rapid relaxation is again followed by an algebraic convergence to the final profile. In particular, in the 
vicinity of y = and y = it (stationary streamlines), it relaxes as slowly as, in this case. The vorticity is 

always zero at the stationary streamlines, the profile in the vicinity is linear (see figure ITS)) . not flat as in the 
case of even perturbations. However, as time goes on, this linear region becomes narrower and narrower, and 
instead the region with the profile proportional to y^ extends towards the stationary streamlines. This indicates 
that, in the large time limit, the proflle is locally parabolic in the vicinity of the stationary streamlines, as in 
the case of even perturbations. The proflle being odd, we remark that such a parabolic profile means that the 
asymptotic vorticity profile is not twice differentiable at the stationary streamlines. 

The final profile obtained from the Laplace transform tools (equation p4p l and the direct numerical simu- 
lations again show excellent agreement (see figure [T5)) . 

The space-time series of |w5,x(2/,i)| and \vs^y{y,t)\, similarly to the case of A{y) = 1, shows a relaxation 
toward zero, illustrating the asymptotic stability of the Euler equations. 

We investigate the velocity asymptotic behavior more precisely. Figure [T7] shows the time series at several 
positions. As shown, their asymptotic forms are \vsx{y,t)\ ~ t~", with a = 1, and \vsy{y,t)\ ^ t^^, with 
/3 = 2, for almost all values of y. Thus we verify that the direct numerical simulation for the asymptotic 
behavior of the velocity perturbation is in good agreement with the theory. Only in the vicinity of y = and 
TT (stationary streamlines), the exponents are changed to a = 1.5 and /3 — 1.5. We recall that we have no 
theoretical predictions for the asymptotic velocity on these stationary streamlines. 

In the last paragraphs, we have compared the mode lus from the direct numerical simulations of the 2D 
(nonlinear) Euler equations with the predictions of the 2D linearized Euler equations. The agreement between 
both is extremely good. There is indeed no visible differences, even for large times. This is in agreement with the 
theoretical discussion in the section [31 namely, the difference is expected to be of the order e^. To summarize, 
we conclude that the asymptotic behavior of the (nonlinear) Euler equations are very well described by the 
quasi-linear theory discussed in the previous sections. 

Evolution of the base flow profile and asymptotic stability. We now consider the evolution of the 
base fiow profile ilo{y,t) (the x-average vorticity, directly related to the a;-average velocity, please see the 
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discussion preceding equation (|43l) page HM ). We define the difference witli respect to tfie initial profile by 
LUo{y,t) — D,o{y,t) + U'{y). We thus have 



The evolution of the base profile ujo{y,t) is due to the nonlinear effects (see Eq. (j43l) page [TS)) . Figure [TH] 
shows this evolution in direct numerical simulations. This illustrates that the perturbation vorticity converges 
extremely rapidly (on times of order t = 15 which correspond to the linear shear times) toward a fixed perturba- 
tion profile. The asymptotic perturbation profile is of order e^. All these are in agreement with the theoretical 
discussions of section |31 

5 Discussion 

In this paper, we have discussed the asymptotic stability of parallel flows for the 2D Euler equations. Our results 
are valid for any flow that have no modes for the linearized dynamics (neither unstable nor neutral ones) . This 
situation is a generic one, as the example of the Kolmogorov flows (section U]) illustrates. An adaptation of the 
present results to the case where the flow has neutral modes would be easy. Our results are valid for base flow 
profiles U (y) with or without stationary points j/o such that U' (yo) = 0- We have emphasized the case with 
stationary points that has not been studied before. 

For the linearized Euler equation, we have proven that Lundgren's ansatz ^ actually describes the asymp- 
totic vorticity field for large times. The asymptotic vorticity field thus oscillates, for each streamline, at a 
multiple of the streamline frequency. The asymptotic vorticity profile depends both on the initial condition in 
a non trivial way, and on the base fiow. The asymptotic vorticity is always strongly affected by the base flow 
structure, in a non-local way, especially when stationary streamlines exist. It is thus unlikely that a descrip- 
tion based on the local shear give a good quantitative description, except may be in a limit or nearly linear 
shear. We have also shown that this asymptotic profile can be computed directly from the resolvent operator 
of the linearized Euler equation (see equation (1341) and figure llip without performing costly direct numerical 
computation of the Euler equations. 

For the linearized Euler equations, we have also proved that the asymptotic velocity field decays algebraically 
for large times (equation [1]) , with exponents a = I and (3 = 2 for any streamlines that are not stationary 
{If (yo) 7^ 0). On the stationary streamlines, we have no theoretical predictions, but we have found numerically 
that two cases exist a = 1 and (3 = 2 for perturbation vorticity fields having the same periodicity as the base flow 
velocity or a = 3/2 and /3 = 3/2 for perturbation vorticity field having the opposite periodicity with respect to 
the base flow. Without stationary streamlines, these results are the same as the classically expected ones. With 
stationary streamlines, these results were unexpected as the effect of the stationary streamlines in oscillating 
integrals could have been expected to give l/^/i contributions. Such contributions cancel out because of a 
self-consistent vorticity depletion at the stationary streamline. This is a new mechanism of vorticity depletion 
at the stationary streamline that we theoretically predict for the linearized 2D Euler equations and prove to be 
self consistent for the 2D Euler equations, and numerically confirm for the 2D Euler equations in this paper. 

This vorticity depletion mechanism occurs due to the effects of the transverse component of the velocity 
perturbation on the background vorticity gradient. This mechanism is thus absent in cases where the background 
vorticity gradient identically vanishes, or for a beta-plane barotropic fiow when the beta effect exactly balance 
the vorticity gradient, a case studied in several papers [Tni [Uj. We think that this last case is not generic as 
the vorticity depletion mechanism exists as soon as the vorticity gradient is not exactly balanced. 

We use the above results to prove that if the perturbation evolves according to the linearized 2D Euler 
equations, the nonlinear term remains uniformly bounded in time, and actually decays algebraically for large 
times. Based on these results, we argue that for the Euler equations (nonlinear), a quasilinear treatment of 
the nonlinear terms is self consistent. This strongly suggests that such a quasilinear treatment of the nonlinear 
term should be valid. This also suggests that the full nonlinear equation converges towards Lundgren's type 
asymptotics for the perturbation vorticity field and to zero for the asymptotic velocity field, extremely rapidly. 

From these theoretical arguments, we then expect that the velocity of parallel flows without unstable or 
neutral linear modes is asymptotically stable : the velocity converges towards a new parallel flow which is very 
close to the initial one, even in the absence of dissipation. The distance between the initial profile and the 
asymptotic one is of order , where e is the order of magnitude of the initial perturbation. 




(49) 
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Direct numerical simulations of the Euler equations close to the Kolmogorov base flow show an excellent 
agreement with the above theoretical predictions. 




Figure 19: A snapshot of the vorticity field for the 2D Navier-Stokes equations with stochastic forces, in a 
statistically stationary regime. The vorticity field is close to a steady state of the Euler equation (here a dipole). 
The fluctuations close to this state are the visible filaments on the figure. One clearly see that such filaments 
are present in between the two vortices, but are absent in the core of the vortices. This is due to the vorticity 
depletion mechanism at the core of the vortices, the points where the angular velocity of the vortices have local 
extrema. 

The theoretical study performed in this paper, could be easily generalized to the study of the asymptotic 
behavior and stability of jets in the context of both barotropic flows in the beta-plane approximations, or 
two dimensional axisymmetric vortices. Many recent works have considered perturbations to two dimensional 
vortices [67l [3l [6l [56l [58] . As far as asymptotic behavior is concerned, following this paper approach, in the 
case of vortices, we argue that a similar perturbation vorticity depletion should occur at any stationary point 
of the angular velocity of the vortex U (r) /r. For instance in the case of a monotonic angular velocity, such a 
stationary point of the angular velocity is located at the core of the vortex r = 0, and the vorticity depletion 
occurs at the center of the vortex. This phenomena has indeed been observed by Bassom and Gilbert [7] (see 
their discussion, and the comment of their flgure 2(b) and 4(a)). They have stated that "We at present lack a 
simple physical explanation of this process whereby vorticity is more highly suppressed than a passive scalar, 
and do not know whether it has applicability beyond the Gaussian vortex". The type of arguments developed 
in the present paper, based on the Laplace tools, suggests that such a vorticity depletion is indeed a generic 
phenomenon, valid for any parallel flow (resp. circular vortex), at the stationary points of the velocity profile 
(resp. angular velocity profile). Mathematically this is due to the regularization of the critical layer singularities 
at the edge of the continuous spectrum. 

This vorticity depletion mechanism also impacts turbulent flows where the perturbations are locally governed 
by the linearized equations. Fig [TOl shows a snapshot of the vorticity field in the dynamics of the stochastic 2D 
Navier-Stokes equations \XU- clearly observes a depletion of the vorticity fluctuations at the core of the 

vortices. This effect in a stochastically forced equation is related to the depletion mechanism in a deterministic 
equation, described in this paper. 

We have treated the linearized dynamics and the asymptotic stability for the case of parallel flows, for the 
2D Euler dynamics. The generalization of such results for more complex cases, for instance flow with separatrix 
and stationary points would also be extremely interesting. The problem is then much more difficult from a 
theoretical point of view, but could be addressed numerically. Also the time dependent situation, by contrast 
to the case of perturbation of steady base flows, is of a large interest. It has been shown numerically that 
interactions with large scales dominate the small scale dynamics ^51139) . in the spirit of rapid distortion theory 
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or quasilinear approaches. This has deep impact on the statistics of the associated turbulence [47]. 

Because both are transport equations by a non-divergent field, there is a very well known analogy between 
the Vlasov equations and the 2D Euler equations. An even closer relation between the 2D-Euler and Vlasov 
equation can also be found: the 2D -Euler dynamics of a localized perturbation (vorticity defect) is actually 
described by a Vlasov equation at leading order [1]. Both the 2D- Euler and Vlasov equations have very similar 
behaviors, including for instance relaxations without dissipation (Orr mechanism or Landau damping) and the 
associated asymptotic stability. The theory of the asymptotic stability of 2D Euler equations is thus deeply 
related to the asymptotic stability of Vlasov equations. We note very recent mathematical results on nonlinear 
Landau damping [46J, a subject related to the asymptotic stability of Vlasov equations. A natural issue, is to 
know if such recent mathematical results [46J, could be generalized to the case of the 2D Euler equations, in 
relation with the results obtained in the present work. 

It is important to notice that a rigorous mathematical treatment of the Orr mechanism for the Euler 
equations, in the spirit of the work by Mouhot and Villani [46J for the Vlasov equation, does not exist yet, even 
in the simplest case of base flow profile U without stationary point. This is a very interesting mathematical 
problem and we hope that some new results may follow the recent proof in the Vlasov case. 

The Laplace tools we have used along this paper is suited for analytical initial data only. An interesting 
question is to know if similar results may hold for non analytic data. We note that for the Vlasov equations, 
analyticity of initial data leads to exponential decay of the perturbation ; but there are counterexamples by 
Glassey and Schaeffer [Ml ISS] showing that there is in general no exponential decay for the linearized Vlasov 
Poisson equation without analyticity, or without confinement. This seems to indicate that analyticity is essential 
for observing Landau damping for the Vlasov equation. We refer to [1^ for a further discussion of this point. 
We guess that similarly, analyticity is an essential ingredient for the Orr mechanism in the context of the 2D 
Euler equations. 

The current work has interesting consequences for the understanding of the kinetic theory of the point 
vortex model in one hand and for the linearized 2D Euler and 2D Navier Stokes equations with stochastic 
forces, when the large scale structures dominate the dynamics on the other hand. These applications were one 
of our motivation for studying the asymptotic behavior of the linearized Euler equation and asymptotic stability 
of the Euler equations. These consequences will be developed in forthcoming works. 
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A) Oscillating integrals 



A-1) General classical results 

Let us consider the asymptotic behavior, for large t, of the integral 

I{t)^ I dxg{x)exp{if{x)t) (50) 

J a 

1. First, we consider the case when / has no singular point. We suppose that / and g are twice differentiable, 
that / has no singular point, i.e. f'{x) ^ for any x, and that either g{a) 7^ or g{b) 7^ 0. Then 



it 



g{h) exp (^t/(6)) _ g(a) exp(zt/(a)) 
/'(&) /'(«) 



(51) 



Namely, I{t) ^ t ^ . This can be easily proved by integrations by part. If g{a) — g{h) = 0, if 5 and / are 
sufficiently differentiable, then the asymptotic behavior can be obtained by further integrations by parts. 

Next, we consider the case when / has a single stationary point. We suppose that / is a real function with 
a single stationary point xq {/'{xq) — 0), that / and g are smooth and that /" (xq) 7^ 0. The asymptotic 
behavior of the oscillating integral (|50p is then given by the classical stationary phase results |27l lU]. 



where e is the sign of /" (.tq). Namely, I{t) ^ t^^/^. 
A-2) Oscillating integrals and the velocity asymptotic expansion 

We apply the general methods of oscillating integrals to the particular case of the computation of the velocity 
from an oscillating vorticity field, like for instance the case given by equation (j38p . We first evaluate the long 
time asymptotics of 

v(y, t)^ f dy' Gfc(y, y')h{y') exp{-tkU (y') t), (53) 



where the function h is assumed to be twice differentiable and where is the Green function for the com- 
putation of the velocity 'v{y) exp {ikx) from a vorticity field uj (y) exp (ikx). We treat explicitly the case of a 
channel —L < y < L. The results are also valid for periodic boundary conditions for y, just by dropping all the 
contributions from the boundary in the expressions bellow. 

Before going into more general discussions, we note that for the case of a linear base flow U (y) = sy, and 
when h{y) is either constant, sinusoidal, or hyperbolic sinusoidal, then an explicit expression for Wk{y,t) is 
obtained (see [19J). 

Clearly, equation (j53p is an oscillating integral. In order to study its asymptotic expansion, we use the 
results of section A-1 of this appendix. 

First, we treat the case of a base flow without stationary point. We use the fact that GA;(y, .) is smooth 
everywhere, except for y' — y (see appendix B). Then we can use the results on the asymptotic behaviors of 
oscillating integrals (section A-1) for both intervals —L < y' < y and y < y' < L independently. Moreover, we 
assume that h (y) is at least twice differentiable. We then obtain 

v^y, t) - - ^ [Gk,x {y, L) h{L) exp {-ikU [L) t) - Gu.x {y, ~L) h{-L) exp {ikU (-L) t) + h{y) exp {-ikU (y) t)] 

t-»l/ks ikst 

(54) 

the first two terms are contributions from the boundaries and the third term is due to the discontinuity of 
Gk,x (2/1 y') for y — y' (see (j59p in appendix B). Here we define s as the minimum shear rate s = min([/' (y)), 
and I is a typical length scale that characterizes the variations of h (y) . Similarly we obtain 



Vy{y,t) 



t»l/ks k^s'^t'^ 



(y, L) h (L) exp {-ikU (L) t) - (y, -L) h (-L) exp {ikU (L) t) - ikh (y) exp {-ikU (y) t) 

(55) 
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We note that there is no contribution of order \/t in this case, because Gk,y{y, y') has a discontinuity of its first 
derivative only, for y' = y (see (j60p in appendix B ; moreover in the case of a bounded domain Gk vanishes at 
the boundaries {Gk,y{y, L) — 0)). 

Next, we treat the case of a base flow with stationary points y„i such that U' {ym) — 0. We assume that 
each stationary point y„i is not degenerated {f"{ym) 7^ 0). We perform the computation for only one of these, 
denoted ?/o, without loss of generality. 

Recall that the classical results on oscillating integrals ([5^ assume the function g (see Eq. [5(1)) to be smooth. 
This is not always the case for us. Indeed, the Green function Gk{y,y') is not smooth for y = y' . However, if 
2/ 7^ 2/0 1 this discontinuity for Gk or for its derivative can easily be handled by dividing the integration interval 
into two subintervals, as has been done in the previous paragraph. Then we conclude that the leading order of 
the asymptotic behavior is still dominated by the contribution of the stationary points. Thus the result (1521) is 
still valid. 

Then, from using ([5^ . we obtain, for y ^ yp, 



v(y,t)^^^G.(y,,o)Myo)y^^^,7(^exp(^— j , (56) 

where is the sign of —kU" (j/o)- We note that the asymptotic expansion has a discontinuity for y = yo, due to 
the discontinuity of the Green function. Through a straightforward generalization of the classical results (|52p to 
oscillating integrals with discontinuous functions g, we can conclude that this discontinuity is regularized over 
a length scale I = y^l/{\kU"{yo) \ t), that decreases with time. 

B) Green functions 

Let us establish the expression and some properties for the velocity Green function G^. The results on the 
discontinuity of the Green function, equations ([55)) and ([50)) bellow, are necessary for the discussion of appendix 
A. We here consider the case of periodic boundary conditions {y 27r— periodic), though the case of a channel 
geometry —L < y < L can be treated similarly, the resultant equations ()59p and ()60p remaining unchanged. 

Let us denote Hk the Green function for the stream function in the x-Fourier space. The x-Fourier transforms 
of w = Alp gives ujk = d'^ipk/dy'^ — k'^ipk- Hk{y,y') is thus solution of 



dy^ 



-eHk^5{y-y') (57) 



where Hk{.,y') = is a 27r— periodic function for any y. 

We note that Hk {y, y'), considered as function of y, is smooth everywhere except for y = y' . For y — y', Hk 
is continuous and has a jump unity for its first derivative: 

dHk , dHk , /-ON 

where F {y'^,y') is the limit of F{y,y') for y' going to y with the condition y > y'. 

Because of the translational invariance in a doubly periodic domain, the set of equations ()57ll58p and the 
associate boundary conditions are invariant under translation. Then clearly, Hk iy,y') depends only on ||y — 
where 

llyll = min |27rn — y\ . 

integer n 

Besides these general properties, an explicit expression to Hk can be found from ()57l) and 

cosh(fc||j/ — 7r||) 



Hk{y,y') = - 



2 sinh(fc7r) 



Using the fact that Hk {y,y') depends only on \\y — y'\\ and ([55)) . it is easily verified that Hk {y,y'), considered 
as a function of y', is differentiable and has a discontinuity in its derivative for y' = y: 

dHk / dHk , . 

-Q^iy^y )-^iy^y ) = i 
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Using V = V A (V'Sz), we have for the x-Foxirier transforms : Wk^x = —di>k/dy and Vfc,^ = ikipk- Thus 

= ( — 5^' 

Then, using the properties of Hk, we note that Gfc is smooth everywhere except for y = y', and that its 

derivative has a jump for y = y' : 

Gk{y,y+)-Gk{y,y-) = {-l,0) (59) 

and 

^{y,y^)-^{y,y-)=ik (eo) 
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Figure 14: The space-time series of the modulus of the perturbation vorticity, \ujs{y,t)\, with initial perturbation 
profile A{y) — siny, and aspect ratio 6 = 1.1. 
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Figure 15: The profiles of the perturbation vorticity modulus, |a;5(?/,t)|, at several times, for the initial pertur- 
bation profile A{y) — siny and aspect ratio 6 = 1.1. As the time goes on, the linear region becomes narrower and 
narrower, and instead the region with the profile proportional to extends towards the stationary streamlines, 
leading to a parabolic profile in the large time limit. 
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Figure 16: The final profile of the perturbation vorticity modulus, |a;5oo(y)|, obtained from both direct numerical 
simulation and theory. The initial perturbation profile is Aijj) — siny and the aspect ratio is (5 = 1.1. The two 
profiles show excellent agreement. 
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Figure 17: The time series of the perturbation velocity components, \vs,x{y,t)\ (a) and \vs^y{y,t)\ (b) at three 
locations, y = (vicinity of the stationary streamHne) (red), y = tt/A (green), and y = tt/2 (blue), for the initial 
perturbation profile A{y) — siny and aspect ratio 6 = 1.1. We observe the asymptotic forms \vs,x{y,t)\ ^ i"", 
with a — 1, and \vs,y{y,t)\ ^ t~^, with /3 = 2, in the almost all the region, in accordance with the theory. Only 
in the vicinity of y = and tt, we observe the exponents a — 1.5 and /3 — 1.5, for which we have no theoretical 
predictions. 41 




Figure 18: The space-time series of the a:-averaged perturbation vorticity, ujQ{y,t). The initial perturbation 
profile is A{y) = 1 and the aspect ratio is i5 = 1.1. 
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